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ABSTRACT 

The apparent presence of large core radii in Low Surface Brightness galaxies has been claimed as 
evidence in favor of warm dark matter. Here we show that WDM halos do not have cores that are 
large fractions of the halo size: typically, r corc /r2oo ^$ 10~ 3 . This suggests an astrophysical origin for 
the large cores observed in these galaxies, as has been argued by other authors. 



1. INTRODUCTION 

The highly successful cold dark matter (CDM) model 
idealizes the thermal motions of dark matter particles 
as negligible on all scales at high redshift. In this 
model, perturbation modes on all scales are gravitation- 
ally unstable, leading to hierarchical structure forma- 
tion in which nonlinear structures such as halos assem- 
ble through numerous mergers. Numerical simulations 
of structure formation within CDM models indicate that 
halos are predicted to have steep central density profiles, 
with logarithmic slopes d log p/d log r ~ —1 on the small- 



est resolved scales (see |Navarro et ah (2010) for a recent 
example) . 

In many DM models, however, the DM temperature 
is nonzero, which can affect the properties of DM ha- 
los in multiple ways. For example, a finite DM tem- 
perature suppresses the abundance of low-mass halos. 
This occurs because, following freezeout of dark matter 
interactions, DM particles freely stream over some dis- 
tance determined by their thermal velocities, and density 
fluctuations on scales below this free-streaming scale n s 
are highly suppressed. Roughly speaking, the smallest 
halos that arise are expected to have masses of order 
M fs = (47r/3)p TO r|, although N-body simulations have 
not definitively ruled out the formation of at least some 
halos below M f ? through non-hierarc hical processes like 
fragmentation ( |Wang fe White||2007| ). 

Besides this suppression of the abundance of halos, 
a nonzero DM temperature also can affect the central 
density profiles of the halos that do form. One elegant 
way to see this is to note that the DM phase space den- 
sity is finite if the DM temperature is finite, and since 
DM is taken to be collisionless, then Liouville's theorem 
guarantees that the phase space density cannot increase. 
Hence, the phase space density is bounded within DM 
halos (Tremaine & Gunn 1979), which implies that the 
central slope ot the D M density profile must vanish (e.g. 



Tremaine et al!p994 1. In other words, halos are expected 
to have central cores if the DM is not cold. Note that 
this effect is caused by the finite DM temperature at the 
time of formation of the halo, and is not due to the trun- 
cation of the small-scale linear power spectrum (Wang & 
|White|[2009l ). 

Therefore, increasing the DM temperature has the ef- 



fect of suppressing the number of low-mass halos, and 
of producing central cores in DM halos. Observation- 
ally, there may be evidence for both of these effects. 
The observed number of Local Group satellite galaxies 
falls well below the thousands of DM subhalos found in 



CDM simulations of halos like our Galaxy's (see Kravtsov 
(20101 for a recent review). In addition, the 21 cm ro- 
tation curves of certain galaxies, in particular low sur- 
face brightness (LSB) galaxies, appear b etter fit by cored 
DM profiles than cuspy DM profiles (see Kuzio de Naray 
et al. (20101 for a recent discussion). For these reasons, 



there has been considerable interest in the literature in 
investigating structure formation in models where DM 
is not perfectly cold. In so-called Warm Dark Matter 
(WDM) models, the DM temperature is chosen to make 
the free-streamin g scale correspond to subgalactic scales, 
r fs ~ Q.l/t- x Mpc ( |Bode et al.||2001[ ). 

WDM models have become increasingly disfavored in 
recent years, in large part because of constraints on the 
matter power spectrum d erived from the Ly man-q for- 



est flux power spectrum (|Seljak et al. 2006). The ob- 
servational support for WDM models from dwarfs and 
LSB galaxies has also eroded, as faint Local Group satel 
lites have been discovered in increasing numbers (e. g 
lokurov et aL|2007| [Toilerud et al.||2008| |Kravtsov| 
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de Naray et al. (2010) have argued that the large cores 
apparently observed in certain LSB galaxies cannot all 
be due to WDM, since the implied central phase space 
densities in these systems are not universal, but instead 
show large variations from object to object. This sug- 
gests an astrophysical origin for claimed detections of 
central cores, and m echanisms to produce s uch cores have 
been proposed (e.g. Governato et al.|[2009 ). 

In this paper, we present yet another argument against 
WDM as the or igin of large cores in halo de nsity profiles. 
As discussed by Kuzio de Naray et al. ( 2010[ ), the inferred 
core radii in several LSB galaxies are large fractions of the 
halo virial radii, r corc ~ 5%r2oo- As mentioned above, 
such large cores do not arise in CDM models, and so we 
might naturally consider producing large cores by raising 
the DM temperature. Making the DM warm, however, 
has the side effect of wiping out small halos, and so it is 
not obvious that raising the DM temperature can gener- 
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ate halos with large r core /r2oo- 

We can, however, use a simple argument to make an 
order-of-magnitude estimate of this ratio. Consider a 
WDM particle of mass m and typical momentum p. Fol- 
lowing freezeout, its momentum redshifts as p oc a -1 , 
so let us write p = mv^/a, where Vg is the veloc- 
ity today at z = (since WDM must be nonrelativis- 
tic today, vo <C 1). The particle's velocity is then 
v(a) = vq/\/vq + a 2 (using units where c = 1). Ne- 
glecting accelerations caused by gravitational potential 
fluctuations, the particle freely streams a comoving dis- 
tance 



vdt v . / a oq 
1/2 — lo S 



(1) 



where we assume vq <C a cq = Q, r /Q. m . This distance 
encloses mass Mf s = (47r/3)p m rf s , and as noted above, 
the smallest halos that form will have masses of order 
Mf s . The virial velocities of these halos at their formation 
epoch a c are 
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H r {s , (2) 



where A 2 oo ~ 200 is the virial overdensity. Now, because 
the thermal velocity u t h ~ ^o /o. c at expansion factor a c is 
nonzero, infalling particles will not fall directly towards 
the halo center, but instead have a nonzero impact pa- 
rameter, and the typical impact parameter determines 
the core size. Naively, we might expect that at forma- 
tion (a = a c ), 
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(3) 



Following formation, any subsequent growth in halo mass 
can only decrease the core radius, while the virial radius 
can only increase. Indeed, even if the mass distribution 
around the halo is static, with no accretion following for- 
mation, the virial radius will still grow in time propor- 
tionally to the expansion factor a, because the expansion 
of the universe dilutes the background mean matter den- 
sity. This is the reason why halo concentrations corre- 
late with halo formation times in CDM c osmologies (e.g. 
Wechsler et al.|2002||Zhao et al.|2003a|b| . Therefore, fol- 
lowing the formation epoch a c , the ratio r coro /r v i r must 
diminish in time at least as fast as a c /a; any mass accre- 
tion will only decrease this ratio even faster. Assuming 
no growth, then at the present time (a = 1) 
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This ratio is maximized by delaying halo formation as 
late as possible. For typical parameters, and setting a c = 
1, this gives r corc ps 10~ 3 r2oo observed today. 

From this simple order-of-magnitude estimate, it ap- 
pears unlikely that WDM models can produce sufficiently 
large core radii to explain LSB galaxies. This argu- 
ment is only approximate, however. To make further 
progress, we have performed calculations of halo forma- 
tion in WDM models. Our results indicate that WDM 



halo cores are broadly consistent with (though typically 
smaller than) the above estimate, which precludes WDM 
as an explanation for the large cores that are claimed to 
exist in certain LSB galaxies. 

2. NUMERICAL METHOD 

In this section we describe our numerical method to 
solve for the self-consistent halo density profile following 
collapse. We eschew conventional N-body simulations, 
since for feasible simulation parameters the core radii will 
typically be unreso lved or (at best) marginally resolved 
(Colin et al. 20081. Since we are interested in studying 
the behavior on scales smaller than the typical resolu- 
tion limits of conventional N-body simulations, we have 
instead employed an alternative approach similar to that 
used by |Lithwick fc Dalal| ( |2010j > . 

We calculate the collapse of isolated peaks in an ex- 
panding universe. To isolate the effects of the nonzero 
WDM temperature during halo collapse, we focus on the 
case of spherical collapse. This problem has been inves- 
tigated previously in the literature, and it is straight- 
forward to show that cold, spherical collapse gives halos 
with ce ntral density profiles behaving as p oc r~ 2 , or 
steeper (|Fillmore fc G oldreich 1984; Bertschinger 1985 
Lithwick & Dalai"'" 1 "' " ' 1 : 1 



2010[ ). Because cold spherical collapse 



2 profiles are 
By Newton's 



is well-understood, any departures from r 
clearly due to the effects of warm collapse, 
theorem, our calculations neglect the effects of the local 
environment of peaks, which cause peaks to collapse non- 
spherically. We know from previous cosmological WDM 
simulations, however, that accounting for the effects of 
local e nvironment does not lead to large core s in WDM 
halos dCoh'n et al |[2008| |Wang fc Wliite|[2009| . 

For a given potential <P{r,t), assumed to be spheri- 
cally symmetric, we integrate the equations of motion 
to solve for the orbit R(t) for each particle. Given the 
orbit R(t), we compute the mass profile deposited by 
each particle. We compute orbits for many particles, 
and sum over all their deposited mass profiles to obtain 
the total mass profile M tota \(r,t), and the total density 
p = (dM t otai/dr)/(4:irr 2 ). Then, we repeat this proce- 
dure, using the newly obtained mass profile, and iterate 
to convergence. 

We initialize this procedure using linear perturbation 
theory. We start with a linear density profile 5(r^) de- 
scribing the initial peak at the starting epoch eiinit — 
(1 + ^init) 1 - Here, tl is a comoving Lagrangian radius, 
to be distinguished from the proper Eulerian radius r at 
subsequent times. We choose the initial peak profile to 
be proportional to the (linear theory) matter correlation 
function, 5(r^) oc £(7~l)- This corresponds to the average 
profile of high peaks in the v — > oo limit ( Bardeen et al. 
1986), and so this profile should be typical of the first 
halos to form in WDM cosmologies. The matter corre- 
lation function depends o n the WDM transfer function, 
which [Bode et al.] (2001) found to be well described by 
the parametnzation 



WDM 



(k) = [l + (akfT 5/u , 



(5) 



where i/k1.2 and a is a characteristic length scale. We 
parametrize WDM models by their free-streaming scale 
rf s , or equivalently the enclosed mass Mf s , so we require a 
translation between rf s and a. We determine the equiv- 
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ale nt free-s treaming length for a given a by matching 
the Bode et al.| transfer function to top-hat smoothing, 
which is given by 



3 

[kRf 



W Tn(k) = 77 ^ w [sin(fci?) - kRcos(kR)} (6) 



for smoothing scale i?, which we take to be equal to rf s . 
We have found that a 5rf s provides a reasonable match 
to the two functions. 

Given a desired free-streaming mass Mf s and halo for- 
mation redshift Z{ orm , we set the initial peak profile shape 
S(r = ainit tl) to be proportional to the correlation func- 
tion (using the appropriate rf s ), and normalize the peak 
height at the starting redshift so that the average inte- 
rior overdensity S = 3r~ 3 L r 5 dr, linearly evolved to 
redshift Zf orm and evaluated at the free-streaming scale 



is 5(r L 



lorm 



Zi orm ) 



1.686, in accordance 



with the spherical collapse model (Gunn & Gott 1972). 
We assume that this linear profile evolves at early times 
according to linear perturbation theory: 



8(r = ar L ,a) = D(a)S(r L ) 



(7) 



where D(a) is the linear growth factor, which for ACDM 
universes with no pressure perturbations may be ex- 
pressed as 



normalized so that D(a = 1) = 1 (] Peebles 1980) 



Note that this procedure is not entirely self-consistent, 
since our use of CDM growth factors neglects the scale- 
dependence in WDM growth factors caused by residual 
free-streaming at late times. Our neglect of this residual 
free streaming when normalizing the initial peak height 
means that our peaks do not actually reach S = 5 C at 
redshift Z{ orm , leading to slight errors in the formation 
epoch. This does not appear to affect our results signifi- 
cantly. 

This procedure specifies the initial overdensity pertur- 
bation at the starting redshift of the simulation. We also 
require the initial velocities for all the particles. These 
velocities have three contributions: the Hubble velocity, 
the bulk peculiar velocity, and a random thermal veloc- 
ity. The Hubble term is of course just given by vh = H r. 
The bulk peculiar velocity may be computed from the 
density profile, using the linearized continuity equation 



6 + V • v = 0, 



(9) 



along with the assumption of potential flow at early times 
(i.e. Vx« = 0). This gives 



v r (r, a) ~ -^rS(r, a)H(a), 



(10) 



where again 5 = 3r~ 3 JT r 2 S dr. In addition to this bulk 
peculiar velocity, for each particle we add a thermal ve- 
locity, drawn from a Fermi-Dirac distribution function 
for WDM temperature T: 



f(p,T)d 3 p 



1 



p 2 dp 



N Q (T) eJ>/ k BT + 1 
where the normalization is given by 

r °° p 2 dp 



N (T) 



o eP/ T + 1 ' 



(11) 



(12) 
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Fig. 1. — Adiabatic orbital evolution. The solid red curve shows 
an example orbit r(t) for a particle from one of our collapse calcu- 
lations. For comparison, the dashed blue curve shows the expected 
behavior for the apoapse under the assumption that the orbit re- 
sponds adiabatically to the deepening gravitational potential. 

We typically begin at redshift z — 100. We sample 6000 
initial radii spaced uniformly in volume, up to a max- 
imum radius chosen to enclose the virialized region at 
z = 0. For each radius, we sample the momentum dis- 
tribution with 500 points and the angular distribution 
with 200 points, weighting the particles according to the 
fraction of initial volume, solid angle, and momentum 
distribution that they represent. 

Given the initial conditions for each particle, we then 
integrate forward the equations of motion using a fourth 
order, variable timestep Runge-Kutta integrator. The 
equations of motion are given by the usual Newtonian 
dynamics: 

d 2 R L 2 , . 

w-w = (13) 

where the angular momentum L = R x v is conserved 
because of the assumed spherical symmetry. 

Given an orbit R(t), the enclosed mass profile de- 
posited by each particle is 



M(r,t) = m p Q[r 



R(t)], 



(14) 



where m p is the mass represented by the particle, and 
0(x) is the step function. Summing over all particles 
gives the total mass profile M total (r, t), and the density 
p{r,t) and potential <f>(r, t) follow easily. Having com- 
puted the mass profile M(r,t) for a given iteration, we 
then use that mass profile in the equations of motion for 
the subsequent iteration. In practice, we bin the mass 
profile using a grid with 350 bins spaced uniformly in 
expansion factor a and 500 logarithmically spaced bins 
in radius, and then linearly interpolate from this grid 
to estimate the mass M at arbitrary times and radii as 
needed for the orbit integrations. 

To expedite this calculation, we have made use of a 
simplifying approximation. For particles deep within the 
halo, whose orbital times are small compared to the Hub- 
ble time, we stop the orbital integration after the frac- 
tional change in the product i? apo x M(i? apo ) over one 
orbit is less than 10 -3 . Thereafter, we assume that the 
orbit evolves adiabatically. Specifically, we assume that 
the radial action J r = § v r dR cx [R a .p M(R ajpo )} 1 / 2 is an 
adiabatic invariant. Given the time evolution of the mass 
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Fig. 2. — Convergence of the density profile. The halo density 
profile at z = 2.7 after iterations f,2,3,6 and 7 are shown in the 
colored curves, while the black curve shows the converged profile 
(12 iterations). 

profile M(r,t), we then easily determine the time evolu- 
tion of the orbital apoapse R apo - As Fig. [I| illustrates, 
adiabaticity is an excellent approximation for these or- 
bits. Similarly, we assume that the ratio of periapse to 
apoapse, i? per i/ Rupcn is also conserved because of con- 
servation of angular momentum. Given i? per i and R apo , 
we then assume that the mass profile deposited by this 
particle is 



M{r, t) 




if r > .R a p 

if i?p 0ri < r < R apo 

if r < i? Der i 



(15) 



which is a good approximation to Eqn. ( 14 ) , except for 
radii very near i? pcr i or i? apo . 

This iterative procedure rapidly converges to a self- 
consistent collapse solution; Figure [2] illustrates one typ- 
ical example. As the figure shows, the interior density 
profile quickly settles into roughly r~ 2 behavior, as ex- 
pected, although there are features at both large radii 
and small radii. The spikes at large radii are caustics, 



a well-known feature of cold spherical collaps e (Fillmore 



fc Goldreich||1984| [Vogelsberger et al.||2009[ ). At small 
radii, discreteness effects of the finite number of particles 
leads to noise in the determined profile. This noise in the 
mass profile enters the equations of motion for particles, 
which affects particle orbits and leads to rapidly devel- 
oping instabilities in the phase-space structure (|Henon | 



1973 Barnes et al. 1986 Henriksen & Widrow 1997). 



These instabilities significantly distort the shape of the 
radial caustics at times following collapse. We have 
checked that if we suppress these instabilities by artifi- 
cially smoothing the potential, the caustics match the 
expected form. Because these instabilities are physi- 
cal, rather than numerical in origin, we have opted not 
to suppress them. Accordingly, our density profiles at 
late times, long after collapse, do not show the expected 
prominent caustics. 

3. RESULTS 

In this section we present results of our calculations. 
In the first subsection, we illustrate the behavior found 
in one typical simulation, and in the following subsection 
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Fig. 3. — Collapsed profile as a function of time. The top panel 
shows the halo density profile p(r), in units of the mean matter 
density p, at various redshifts before, during, and after collapse at 
z ~ 2.7. The middle panel shows the 3-D velocity dispersion. The 
bottom panel shows the pseudo-phase-space density p/o" 3 . 
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Fig. 4. — Radial behavior. Depicted are the radii r and radial 
velocities v r for a subset of particles at z = 2.7. The color for 
each point indicates each particle's turnaround radius. The vertical 
dotted lines indicate r core and r200 at this redshift. 

we describe how the behavior changes as we vary several 
physical parameters. 

3.1. Anatomy of a WDM halo 

In this subsection, we describe the detailed results of 
one of our simulations. The behavior found for this halo 
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is representative of our simulations. For concreteness, we 
use a WDM temperature To/to = 1.3-1Q -8 in units where 
c = 1, which gives a free-streaming scale rf s ~ 100 
kpc, and normalize the initial peak amplitude so that 
halo formation occurs near redshift z — 3. Figures [3] 
and [3] show the results. Figure [3] plots the density profile 
p(r), the 3-D velocity dispersion <r(r), and the pseudo- 
phase-space density p/cr 3 , as a function of time, while 
Fig. [4] illustrates a snapshot in time of the radial dy- 
namics. At early times, prior to collapse, the density 
evolves perturbatively, so that the p(r) profile is similar 
to the original linear density perturbation 8(r), simply 
growing in amplitude. The phase-space density remains 
very homogeneous before collapse. Near the time of col- 
lapse, there is smooth infall towards the halo, and both 
the density and velocity dispersion rise in concert to keep 
pjo nearly constant. Orbits do cross in the infall region, 
however, since particles with different thermal velocities 
fall into the collapsing halo at different rates. 

The initial collapse of the halo produces a roughly r~ 2 
density profile, due to our assumption of spherical sym- 
metry. This breaks to a shallower p ~ const behavior at 
the core radius r corc , where the infalling particles reach 
periapse. Following periapse, the particles splash back 
outwards with positive radial velocity. The splash-back 
surface defines the outermost caustic, where both the 
density and velocity dispersion jump precipitously. The 
density jumps at the caustic due to a pile-up of particles 
with similar apoapses. The velocity dispersion jumps be- 
cause outside the caustic, particles are all falling inwards, 
whereas inside the caustic there is both outwards and in- 
wards motion. Inside the caustic, the velocity dispersion 
remains roughly isothermal, at a value near the halo's 
virial velocity. 



Because p 



and a is nearly constant with radius, 



the pseudo-phase space density shows a nearly power- 
law behavior over much of the virialized region. Outside 
the outermost caustic, pjo 3 remains nearly identical to 
the phase space density of the unperturbed material. At 
the caustic, p/a 3 falls steeply, due to the sudden increase 
in velocity dispersion. Towards smaller radii, p/a 3 rises 
smoothly, close to r~ 2 . It is important to stress, however, 
that pjo 3 is not a good proxy for the actual phase-space 
density over much of the halo's interior. The velocity 
dispersion tensor is highly anisotropic, in the sense that 
radial velocities are much larger than tangential veloci- 
ties, as can be seen from the predominantly radial orbits 
shown in Figures [l] and HI For this reason, er « ay, and 
so a 3 ~ a 3 ^> a r aga,p. Only near the core radius does 
the velocity dispersion become close to isotropic. 
Many of these features are similar to what is found 



in cold spherical collapse (Fillmore & Goldreich 
Bertschinger 1985 1. The most obvious difference be 



1984 



tween 



warm and cold collapse is the presence of a core radius, 
caused by the orbits' inability to reach r = due to their 
nonzero angular momenta. We estimate the core radius 
by fitting the density profile to the functional form 



(16) 



1 



2/a ' 



The parameter a controls how sharply the profile breaks 
from r~ 2 behavior to constant density, and typically our 




Fig. 5. — Time dependence of halo properties. The upper, middle 
and lower panels show r corc , r200, and the ratio r COTe /r 200 respec- 
tively, as a function of redshift. As the halo grows in mass over 
time, the virial radius grows, while the core radius shrinks due to 
adiabatic contraction. The ratio r core /r200 is therefore maximized 
at the time of halo formation, and only diminishes thereafter. 

simulations give a ~ 1 — 2. We are mainly interested, 
however in the core radius. As Figure [4] shows, r corc co- 
incides with the typical location of the periapse of the 
infalling particles. Both the density and pseudo-phase 
space density plateau at r core , the latter saturating at a 
value near the phase space density of unvir ialized mate- 
rial outside the halo as expected from the Tremaine & 
Gunn| (| 1 9 79 1 ) bound. 

Figure [5] illustrates the evolution of the halo structural 
parameters over time. Following formation, the halo con- 
tinues to accrete matter and steadily grows in mass, at a 
rate determined by the initial linear overdensity profile. 
By definition, this growth in M200 means that r2oo grows 
as well. Note, however, that r coie decreases as the halo 
grows, due to adiabatic contraction of the orbits as the 
halo potential deepens over time. Since r corc only shrinks 
in time, while r2oo grows in time, the ratio between the 
two is clearly maximized at the time of formation of the 
halo. 

Lastly, in Figure [6] we plot the breakdown of the con- 
tributions to the mass profile from various radii. The 
left panel decomposes the particles into shells of initial 
radius at z = 100, while the middle panel uses bins of 
turnaround radius r ta . At large radii, r 3> r core , the 
mass is dominated by recently infalling particles that 
originated at large Lagrangian radius and have large 
turnaround radii, similar to the behavior fou nd in cold 
spherical collapse (Fillmore & Goldreich 1984). However, 
this changes on scales of order the core radius. Near 
r CO rej many different shells spanning decades in radius 
contribute comparably to the density. When we bin the 
particles based on initial, Lagrangian radius, we can see 
that each shell has a different core radius, roughly scaling 
as r core cx tl, as may be expected from the arguments 
given in the introduction. When we bin the particles 
based on their turnaround radii, we can see that the typ- 
ical apoapse for each shell scales like the turnaround ra- 
dius. For shells with r apo > r corc , the enclosed mass 
profile behaves as 
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Fig. 6. — Breakdown of the mass distribution. The left panel shows the contribution to the mass density p(r) at z = 2.6 originating 
from shells at various initial radii at z = 100. The middle panel shows a similar breakdown, instead binning particles based on their radii 
at turnaround. The right panel shows the distribution of initial radii (blue dashed) and turnaround radii (red solid) for all particles with 
r < r coro at z = 2.6. 



and since 



oc 



tl for each shell, we see that inside 



the core radius, the mass profile deposited by each shell 
of width dlogrL scales as M s h e ii oc r 3 x (rL/r apo ) oc 
r 3 x (r L /r ta ). (For shells with r apo < r corc , the en- 
closed mass simply behaves as M s h e ii(r) oc (rLr/r corc ) 3 , 
of course.) For cold collapse, there is a one-to-one re- 
lationship between the turnaround radius r ta and the 
initial radius, that depends on the initial linear density 
profile of the peak collapsing to form the halo. Roughly 
speaking, if locally the linear density has slope 7, in 
the sense that S oc r7 7 , then r ta oc ?\ +7 (Lithwick & 



Dalai 2010). The tight relationship between tl and rt a 
degrades somewhat for warm collapse, but we can still 
use the same basic scaling. For shells inside the free- 
streaming scale, tl < rf s , the peak profile is quite flat 
with local slope 7 ~ 0, and so r ta oc r^. For these shells, 
-MshciiM becomes roughly independent of the shell's ini- 
tial radius for r < r COIC , as is seen in Fig. [6] For larger 
radii, tl > rf s , the slope of the initial profile becomes 
nonzero, 7 > 0, implying that shells originating from 
large radius become subdominant inside the core radius. 
Thus, the mass inside the core radius is dominated by 



shells with tl < ff s and r a 



> r c 



3 , and receives roughly 



equal contributions per decade within this range, as is 
seen in the right panel of Figure [6] 

3.2. Dependence on physical parameters 

Having established the basic features of the halo struc- 
ture, we now explore the physics that sets those proper- 
ties. The two main differences between WDM and CDM 
are the cutoff in the power spectrum, and the relic ther- 
mal velocities of DM particles at the time of halo forma- 
tion. Both of these differences influence the size of halo 
core radii. To disentangle the different effects, we have 
performed collapse simulations in which we hold fixed 
the linear density profile of the initial peak, but vary the 
WDM temperature. This corresponds to holding fixed 
the halo assembly history, but varying the random mo- 
tions near the time of collapse. The argument given in 
section 1 (e.g., Eqn. (|3|) would predict that r coro would 
scale linearly with temperature, and our calculations ap- 
pear consistent with this, as shown in Fig. [7] As we vary 
the temperature, the overall assembly history and struc- 
ture of the halo remains unchanged (e.g., the location and 
height of the caustics), however the core radius varies. 
We find that a simple linear scaling, r corc oc T, appears 
consistent with our simulations. We note, however, that 
this linear behavior breaks down at very high temper- 



atures, when the particles' random velocities become of 
order the Hubble velocity at the time of halo collapse. In 
this regime, the thermal motions are no longer a small 
perturbation to the particle dynamics, and the overall 
collapse of the halo is significantly modified, unsurpris- 
ingly. Of course, such calculations are not self-consistent: 
the large random motions that modify halo collapse at 
low redshift would have erased the initial linear density 
perturbations responsible for the halo, at a higher red- 
shift. 

This establishes that at fixed initial peak profile (i.e. 
fixed halo assembly history), the core radius scales lin- 
early with DM temperature. It is inconsistent, however, 
to hold fixed the initial peak profile while the tempera- 
ture is varied, since the random thermal motions of DM 
particles erase structure at high redshift and modify the 
peak profiles in the linear regime of structure formation. 
Therefore, we next explore how the core radius behaves 
as we self-consistently vary both the WDM temperature 
and the initial peak profile. We know that rf s scales 
roughly linearly with T, and we have just seen that at 
fixed rf s , the core radius r core also scales linearly with T. 
Therefore, if r corc were independent of the halo assembly 
rate, then both r corc and rf s would scale linearly with 
T, and the ratio 7" corc /r2oo would be independent of the 
WDM temperature, as we argued in Section 1. Figure 
[8] shows that this behavior is not confirmed by our sim- 
ulations, however. The figure shows the results of sim- 
ulations using temperatures 5 times larger, and smaller, 
than our fiducial calculation. The red solid curve shows 
the density profile for n s — 20 /i -1 kpc, orange dashed 
shows our fiducial run with n a = 100 ft -1 kpc, and blue 
dotted shows results for j-f s = 500 ft. -1 kpc. In all three 
cases, we set the initial peak height so that collapse will 
occur near z — 3. As expected, r2oo scales close to lin- 
early with rf s : the three simulations give r2oo = 1.2, 5.5, 
and 22.6 kpc at the formation redshift z = 2.7. However, 
r core does not scale linearly with rf s : the three simula- 
tions give r core =0.6, 4.8 and 43 pc respectively. The ratio 
r C ore/?"200 is n °t independent of n s , but instead behaves 
1/2 

roughly as rJ over the range that we have considered. 
Evidently, the core radius depends not only on the DM 
temperature at the time of halo formation, but also upon 
the halo assembly rate. 

Lastly, we examine the dependence of the core radius 
on the halo formation time. We do so, simply by ad- 
justing the height of the initial peak, holding fixed the 
WDM temperature and the radial shape of the peak pro- 
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Fig. 7. — Temperature dependence of the halo profile for a fixed initial peak profile. (Left) The black solid curve shows the halo profile at 
redshift 2 = 2.7 for T = To, where To is the WDM temperature consistent with the linear overdensity of the initial peak. For comparison, 
the red, orange and blue curves show the profile for T = O.lTo, T = 0.3T) and T = O.6T0 respectively. (Right) The points show the fitted 
core radii for the profiles depicted in the left panel, using Eqn. (16) , while the blue line shows a simple linear scaling r core oc T. 




Fig. 8. — Effect of WDM temperature. Density profiles at 2 = 2.7 
for different initial peak profiles consistent with WDM temperature 
that produce rf a = 20, 100 and 500 kpc/h are plotted in dashed 
red, solid black and dotted blue lines respectively. 

file. Figure [9] shows one example, comparing our fiducial 
simulation (with Zf olm — 2.7) with a run using the same 
WDM temperature, but with initial peak height a fac- 
tor of 4 larger. The later simulation has Zf oxm = 13.6. 
The halo masses of the two simulations are similar, as 
expected: r 2 oo = 1-3 kpc for the z = 13.6 halo, com- 
pared to r2oo = 5.5 kpc for the fiducial z — 2.7 halo. 
The core radius for the earlier-forming halo is r coro =2.8 
pc, compared to 4.8 pc for the fiducial halo, so that the 
ratio r core /r2oo changes by a factor of 2.8. The sim- 
ple argument given in section 1 would have predicted 
that r corc /r2oo would scale as (1 + Zf oml ) 1 '' 2 , whereas 
our simulation appears more consistent with a scaling 
r C ore/?"200 oc (1 + Zf oml ) 2 '' 3 . This is only based on one 
comparison, of course. This is the result at the forma- 
tion time; at z — the ratio r core /r2oo would be smaller 
by at least a factor of 1/(1 + Zf OTm ), as we argued earlier. 

4. DISCUSSION AND SUMMARY 

We have investigated the formation of halos in warm 
dark matter cosmologies. Our study of spherical collapse 
of WDM halos indicates that core radii do indeed arise in 
these cosm ologies, as expected from simple phase-space 
arguments (Tremaine & Gunn 1979). However, we find 
that WDM core radii are genetically small, typically of 




Fig. 9. — Effect of initial peak height. The solid black line shows 
the profile at collapse time 2 = 2.7 for a peak height normaliza- 
tion of 5(rL = 100 /i -1 kpc, z = 3) = 1.686. In contrast, dashed 
blue shows the profile at collapse time 2 = 13.6 for a peak height 
normalization of <5(rL = 100 kpc, 2 = 15) = 1.686. 

order 10 -3 of the halo virial radius at the time of forma- 
tion, and considerably smaller following formation. This 
is for halos forming at the cutoff scale; higher mass ha- 
los will have substantially smaller r coro /r 200- We have 
investigated the dependence of the core radius on vari- 
ous physical parameters such as the WDM temperature, 
halo formation redshift, and halo mass. For the allowed 
range of WDM temperatures (e.g. Mf s < 1O 9 M ), the 
core radii of halos observed at z — are generically ex- 
pected to be far smaller than the core sizes measured in 
certain LSB galaxies, with r core /r2oo ~ 0.05. 

Our calculations have all assumed spherical symme- 
try, whereas halo formation in both CDM and WDM 
cosmologies is highly nonspherical. We would argue, 
however, that our conclusions regarding WDM cores are 
likely to be valid for non-spherical collapse as well. One 
line of evidence supporting this argument is the fact that 
the pseudo-phase-space density profiles of our halos are 
quite similar to the profiles of halos in CDM simulations . 



with p/a 3 ex r 2 , roughly speaking ( Ludlow et al.|2010 ) 



This similarity is presumably a consequence of the virial 



theorem, which ensures that p/a 



For WDM 



halos, the power-law rise of the pseudo-phase-space den- 
sity towards small r saturates when p/a 3 approaches the 
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Tremaine-Gunn bound, and this saturation will occur 
for both spherical and nonspherical collapse. Now, we 
would expect the value of the pseudo-phase-space den- 
sity to be similar at the virial radius for both spherical 
and non-spherical collapse, since the halo mass and virial 
radius are by definition the same in the two cases, and so 
p ~ M/4nr 3 and a ~ (GAf 2 oo/ r 20o) 1 ^ 2 wu l be similar for 
the two cases. We have noted that p/a 3 rises as roughly 
r -2 inside the virial radius in spherical and nonspher- 
ical collapse, and in both cases the core radius occurs 
where p/a 3 approaches the Tremaine-Gunn bound. So 
we have good reason to believe that halo core radii will 
not be significantly larger for nonspherical collapse than 
for spherical collapse, just because the p/u 3 profiles ap- 
pear similar. 

The other possible loophole in our argument is that we 
have assumed that no halos form below the cutoff scale 
in the power spectrum. N-body simulations have not 



yet conclusively determined whether or not halos with 
M <C Mi s arise in WDM cosmologies, due to numerical 
difficulti es associated with sim ulating truncated power 
spectra ( |Wang fc White 2007). We have begun inves- 
tigating this issue, and our preliminary results indicate 
that halos may form below the cutoff scale, though it is 
unclear whether they can form in sufficient numbers to 
account for observed LSB galaxies. 

Our results indicate that warm dark matter cosmolo- 
gies cannot produce halos with core radii large enough to 
account for the density profiles of observed LSB galaxies. 
This would suggest that the origin of these observed cores 
lies within astrophysics, rather than particle physics. 

We thank Manoj Kaplinghat, Mark Vogelsberger and 
Larry Widrow for useful discussions. FVN thanks CITA 
and the CfA for hospitality during the course of this 
work. ND is supported by CITA, and FVN is supported 
by CSIC pre-doctoral grant JAE 2008. 
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